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Abstract 

The energy density of a scattering soliton solution in Ward's integrable chiral model 
is shown to be instantaneously the same as the energy density of a static multi-lump 
solution of the CP"^ sigma model. This explains the quantization of the total energy in 
the Ward model. 
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1 Introduction 



The Ward model (Ward 1988) is one of the most interesting models of sohton dynamics. It is 
defined in two space dimensions and is almost fully Lorentz invariant. The field equation can 
be regarded as the (2+l)-dimensional Bogomolny equation for Yang-Mills-Higgs fields, with 
gauge group SU(2). This formulation is first order in time derivatives and is fully Lorentz 
invariant, but there is no natural Lagrangian and hence Noether's theorem is not available to 
construct conservation laws. 

An alternative formulation, due to Ward, involves fixing the gauge, and this requires a 
choice of spatial direction. It leads to a field equation for an SU(2)-valued "chiral" field J 
which is second order in time derivatives, 

dt{j-'dtJ) - d,{j-'d,j) - dy{j-'dyj) + [J-'dyJ, J-^dtJ] = . (l.l) 

This equation is a variant of that of the standard chiral sigma model in 2+1 dimensions for J, 
where the final commutator term is absent. In this formulation, Lorentz invariance is broken to 
(l+l)-dimensional Lorentz invariance. However, there is a Lagrangian, and hence a conserved 
energy and a conserved momentum in the ^/-direction. One might imagine that by unfixing 
the gauge one could recover a complete conserved energy and momentum vector, but it is not 
clear how this should be done. 

A remarkable feature of the Ward model is that it is integrable, both from the point 
of view of the (2+l)-dimensional Bogomolny equation and from the point of view of the 
equation for the chiral field J. The integrability ultimately arises from the fact that the 
Bogomolny equation is a reduction of the (2+2)-dimensional self-dual Yang-Mills equation, 
which can be treated by twistor methods. Because of the integrability there is a holomorphic 
structure, an infinite set of conservation laws (loannidou & Ward 1995), including the energy 
and momentum mentioned above, and a procedure for obtaining general solutions in terms of 
initial data using spectral methods and via the solution of a Riemann-Hilbert problem (Fokas 
& loannidou 2001). For us, the most useful consequence of integrability is that there are 
families of algebraically explicit solutions representing multi-soliton scattering processes. 

These soliton solutions are constructed from auxiliary functions of a spectral parameter 
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A, with poles. If there are a finite number of simple poles then the solution describes a finite 
number of solitons moving along straight lines in the plane at constant speed, and they do not 
scatter. 

Ward (1995) noticed that if there is a double pole in A then the solution represents scat- 
tering solitons. The double pole can be obtained (with care) in the limit that simple poles 
coalesce. The construction is technically a bit more difficult than the simple pole case. It is 
related to the construction of static solutions of the planar chiral sigma model, in terms of 
so-called unitons (Uhlenbeck 1989). Whereas the simple pole solutions correspond to single 
unitons, the double pole solutions correspond to 2-unitons. The original example of a double 
pole solution was generalised by loannidou (1996), and a number of interesting multi-soliton 
scattering solutions were found. The general solution of this type depends on an arbitrarily 
chosen pair of rational functions f{z) and h{z), where z — x + iy is the complex coordinate in 
the spatial plane. Recently, a very general analysis of multi-pole, multi-soliton solutions has 
been given (Dai & Terng 2004). 

In this paper we shall make some observations concerning the energy density and total 
energy of the general double pole solution discussed by Ward (1995) and loannidou (1996), 
which has no net momentum. Our proofs of these observations are not deep, and we expect 
that a better understanding is possible. 

Our first observation is that the total energy of the double pole solutions is quantized 
in units of Stt. This was actually noted earher by loannidou & Zakrzewski (1998) in a few 
examples, but no proof was given. It is physically quite a surprising result, because the 
solutions represent solitons in motion and part of the energy is kinetic. We are familiar with 
quantized energy for multi-solitons at rest, if they satisfy a Bogomolny equation, but usually 
any kinetic energy raises the total energy above the quantized value (see e.g. Manton & 
Sutcliffe 2004). We shall clarify below the number of energy units that a particular solution 
has - the relation to the number of solitons is not particularly simple. 

Our second observation is that the energy density of the solutions (including both the 
gradient and kinetic energies of the field), calculated at any instant of time t, is the same as 
the energy density of a static CP^ sigma model solution (also known as a CP^ multi-lump). 
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The CP^ multi-lump has t as a parameter, but is not regarded as dynamical and has no kinetic 
energy. We have no deep understanding of this relationship between the SU(2) Ward model 
solutions and the CP^ multi-lumps; we were led to it through the algebraic form of the energy 
density in a number of examples. Our proof that the energy density is the same is based on 
direct calculation. 

The first observation follows from the second. It is well known that the total energy of 
a static CP^ multi-lump solution is quantized (Zakrzewski 1989), and we can determine how 
many units of energy there are in terms of the degrees of h, f, and the z-derivative of /. 

2 Double Pole Solutions of the Ward Model 

The simplest solutions of the Ward model arise from a function of A with a single pole 
where // is a constant, and R{z,z) is a 2 x 2 matrix depending on a rational function f{z). 
For II — i this solution is static and is equivalent to the multi-lump solution of the CP^ sigma 
model given by the rational map f{z). The topological charge of the lump is the degree of 
/, and the energy in the Ward model is the same as in the sigma model, namely SttA^. 

The next simplest solutions in the Ward model are those based on a superposition of simple 
poles in the A-plane. Let us denote the pole locations by /ii, 112, ■ ■ ■ , A*n) which are necessarily 
distinct. Associated with these are rational functions fi{z),f2{z),...,fn{z). The function 
friz) describes a multi-lump which moves at a constant velocity determined by /ir- Since 
the pole locations arc distinct, so arc the velocities, so at most one of the multi-lumps is at 
rest. Although the multi-lumps are in motion, they do not scatter off each other in this type 
of solution. There is a simple formula for the total energy, which varies continuously as the 
parameters /ir vary, so the energy is not quantized. 

Scattering soliton solutions are obtained by considering a limit of the solutions above, in 
which two or more of the pole locations fir are brought into coincidence. The simplest case is 
where n = 2 and the parameters fii and 112 are brought into coincidence at i. If the limit is 
taken appropriately, a double pole at X — i results, and it is in this case that the momentum 
in the y-direction vanishes and the energy is quantized as a multiple of Stt. Solutions with a 
double pole at a different location are Lorentz boosted and have a different energy. 
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It was shown by Ward (1995) that for these double pole solutions the SU(2) matrix J 
has a factorized 2-uniton form. We recall now the expression for J. It depends on two 
rational functions f{z) and h{z) (which arise from the functions fi{z) and f2{z) as the limit 
of coalescing poles is taken). Explicitly 

where / is the 2x2 unit matrix, and qi and q2 are the 2-component row vectors 

gi = (l + |/|2)(i,/)-2z(t/' + /^)(/,-l), (2.2) 

q2 = (1,/). (2.3) 

q2 is just a function of z, but qi and hence J are functions of z, z and t. f denotes ^. Because 
of the time dependence of gi, the solution describes the scattering of solitons. 

loannidou (1996) explored a number of examples of the 2-uniton solution 1)2.11) . The 
simplest exhibiting soliton scattering is the solution with f{z) = z and h{z) = z^ (Ward 
1995). The scattering can be deduced from the zeros of the expression tf' + h = t + z'^, noting 
that these are located at z = ±a/^. The zeros are on the x-axis for t < 0, pass through the 
origin at t = 0, and are on the y-axis for t > 0. A slightly more complicated example is with 
f{z) = z and h{z) = z^ . Here three solitons scatter with equilateral triangular symmetry. 
Solutions with /i = are also possible. If f{z) = z and h = 0, for example, then there is 
circular symmetry for all t, with a ring-like soliton contracting and then expanding. 

It is rather curious that these solutions (and also their energy density, which we investigate 
in the next section) appear to respect the rotational symmetry of the plane, despite our earlier 
remarks about how Lorentz invariance is partially broken. It is possible that these solutions 
have zero momentum and the same energy no matter which direction is chosen to fix the gauge 
of the Bogomolny equation, but we have not established this. 



3 Energy of 2-Uniton Solutions 

The energy density of the field J in the Ward model is 

Sj = -^TT{{j-'dtjf + {J-'d,jf + {r'dyjf) . (3.1) 
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This is actually the same as in the standard Lorentz invariant sigma model - it is the field equa- 
tion of the Ward model that breaks the Lorentz invariance. If we use the complex coordinate 
z — X + iy then 

£j = -^Tr a J-'dtJf + 4iJ-%J)iJ-%J)) . (3.2) 

The contribution —^TY{J~^dtJy is the kinetic energy density; the remaining gradient terms 
can be thought of as the potential energy density. The total energy is £■ = J Sj cPx. 

Our main result is that this energy density is the same as that of a CP^ sigma model lump 
solution. In the CP^ sigma model (Zakrzewski 1989), the basic field is a 4- vector of complex 
functions in the plane 

V ^ {V,,V2,Vs,V^) , (3.3) 

whose energy density and other properties are unaffected by multiplying all components by a 
common function W. The components of V should not simultaneously vanish. For a finite- 
energy multi-lump solution, all components of V must be rational functions of z. V can then 
be converted to a 4- vector of polynomials in z, with no common root, by multiplying through 
by the common denominator. (For an anti-lump, one takes rational functions of z instead. 
There are also higher energy saddle point solutions which depend on z and z, but these do 
not concern us here.) 

The energy density in the sigma model is 

V II' I'' - V V'\^ 

m\^+\v2\'+\vs\^ + wr' ^ ^ ^ 

where V/ = The topological charge N of the lump is the degree of V, degV, which is the 
highest degree among the polynomials comprising V (after clearing denominators), and the 
total energy J Sy d^x is SttA^. This follows from a Bogomolny argument which reduces the 
energy to Stt times the integral over the plane of the pull-back of the 2-form generating the 
integer cohomology ring of the manifold CP^. 

The multi-lump that corresponds to the 2-uniton solution of the Ward model has the 
specific form 

V^{2{tf + h),f,V2f,l). (3.5) 
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Note that t appears here exphcitly as a parameter, but we calculate the energy density at each 
instant as if the multi-lump were static, using the formula ()3.4|) . 

To find the total energy, we need to determine the highest degree of the polynomials that 
occur in ()3.5|) when one clears denominators. Let deg/ and degh denote the degrees of the 
rational functions / and h. Generically, / = ^ and h = ^, where r and s are polynomials of 
degree deg/, u and v are polynomials of degree degh, and s and v have no common roots. 
Then substituting in (j3.5j) and multiplying by the common denominator s'^v, we find 

V = {2t{sr'v — rs'v) + 2us^ , r^f , v^rsf , s^v) , (3.6) 

so = 2 deg / + deg h. However there are plenty of non-generic possibilities. For example, 
if / and h are polynomials, then is the greater of 2 deg / and deg h. Generically, is also 
the number of zeros of tf + h, but again there are non-generic examples where this is not the 
case. 

Now consider the energy density for V of the form ()3.5p . Since time derivatives play no 
role here, let us denote 2{tf' + h) by g. One finds after some algebra that 



(|^?P + (1 + |/|2)2)2 ■ (3.7) 

We now wish to demonstrate that Ej = Ey. Let us consider the kinetic contribution to £,j 
first. Note that J = AB (the product of the two uniton factors in 1)2.11) ). and that A^^ = A 
and B-^ = B. Note also that B is independent of t. Hence Ti{J'^dtJy = -Ti{dtAdtA). 
Further simplification is possible by expressing qi as a linear combination of the everywhere 
orthonormal, time independent row vectors (1 + |/P)"2 (1, /) and (1 + |/P)-2 (/, -1). Their 
coefficients are, respectively, (1 + |/P)^ and — 2i(l + \f\'^)~^g. After a modest calculation 
involving these coefficients, and noting that the time derivative of g is 2/', one finds that 

--Tr(J dtJ) - (|^|2 + (i + |j|2)2)2' (3.8) 

which is identical to the contribution to Sy from the third term in the numerator of ()3.7|) . A 
rather more involved calculation, using MAPLE, leads to the result that the gradient energy 
is 

-2MJ-%J)iJ-%J)^'-l^l±m^Mf^^ (3,9) 



which is identical to the contribution of the first and second terms in the numerator of ()3.7|) . 
In total, we see that £j = £vi and note that in £j we include time derivatives, but in £y we 
do not. 

Let us illustrate these formulae with a pair of examples. The first is the two-soliton solution, 
with f{z) = z and h{z) = z^. The energy density is 

(5r4 + 2r2 + l + 4t2 + 8t(x2-?/2))2 ' 

where = x"^ + y"^. The energy is peaked at two points on the x-axis for t < 0, forms a ring 
around the origin at t = 0, and is peaked at two points on the y-axis for t > 0. See Ward 
(1995) for figures showing the energy distribution at different times. The total energy is IGvr, 
because the highest degree of the polynomials in \^ is 2. This is also easily verified by direct 
integration of the circularly symmetric energy density at t = 0, and at other times follows 
from energy conservation. 

A second example is f{z) = z and h{z) = z^. Here the energy density is 

+ 16r6 + 19r^ + 2r2 + 1 + 4t2(2r2 + 1) - 8t(r2 + 2)a;(a;2 - 3?/2) 

£ = 16 7^ , (3.11) 

(4^6 + ^4 + + 1 + 4t2 + 8ta;(a;2 - 3y2))2 

(which simplifies and corrects formula (21) of loannidou (1996)). This energy density exhibits 
the scattering of three lumps with triangular symmetry, whose locations are approximately at 
the zeros of z^ + t (see loannidou (1996) for figures). The total energy is 247r. 

We have also investigated other cases where / and h are rational functions, but not poly- 
nomials. Figure H presents an example of generic type, where f = 1/z and h = l/{z — 1). 
The expression for the energy density is complicated, but it is easily seen that degl^ = 3, im- 
plying that the total energy is 247r (as confirmed numerically). In addition, the motion of the 
configuration is interesting as seen from the zeroes of tf + h, that is, the zeroes of z"^ — tz + 1. 
These show the locations of a pair of soliton structures whose motion is not so simple. For 
t < and t > A both zeros are real, but for < t < 4 the zeros are a complex conjugate pair. 
At t = and t = 4 there are repeated zeros at 2; = and z = 2, respectively. In Figure Q we 
can see solitons approaching along the x-axis at early times, colliding at t = and scattering 
at right angles, then moving apart on circular orbits and colliding again at t = 4, and finally 
separating along the x-axis at late times. 



An example of the non-generic type (because the denominator of h is the square of the 
denominator of /) is f{z) = 1/z and h{z) = {z — l)/z'^. Here, the energy density is 

+ 2r2 + 3 + SrV^ + Ar^ - Srr^x ,^ ^ ^, 

t = lb o , (3.12) 

(r4 + 6r2 + 1 + 4r2 - 8rx) 

where t = t + 1. This is symmetric under the simultaneous reflections x —>■ —x and r — r. 
Figure 121 presents the evolution of the configuration for times around t = —1. The total energy 

is 1677. 



4 Generalisations 

More complicated multi-pole solutions of the Ward model can be constructed. For a recent, 
general analysis of these, see Dai & Terng (2004). A triple-pole example, for which J has a 
3-uniton form, is presented in loannidou (1996). In this, and a few further examples, we have 
verified that the total energy is an integer multiple of Svr. We suspect that the instantaneous 
energy density can again be identified with that of a static sigma model solution, but now for 
CP^. Clarification of this, and its generalisation, would be desirable. 

Acknowledgement 

N.S.M. thanks the Mathematics Division of the Aristotle University of Thessaloniki for hos- 
pitality. 



8 



References 

Dai, B. & Terng, C.-L. 2004 Backlund transformations, Ward solitons, and unitons. Preprint 
math.DG/0 405363 

Fokas, A. S. & loannidou, T. A. 2001 The inverse spectral theory for the Ward equation and 
for the 2+1 chiral model. Commun. Appl. Anal. 5, 235-246. 

loannidou, T. 1996 Soliton solutions and nontrivial scattering in an integrable chiral model in 
(2+1) dimensions. J. Math. Phys. 37, 3422-3441. 

loannidou, T. & Ward, R. S. 1995 Conserved quantities for integrable chiral equations in 2+1 
dimensions. Phys. Lett. A 208, 209-213. 

loannidou, T. & Zakrzewski, W. J. 1998 Solutions of the modified chiral model in (2+1) 
dimensions. J. Math. Phys. 39, 2693-2701. 

Manton, N. S. & Sutcliffe, P. M. 2004 Topological solitons. Cambridge University Press. 
Uhlenbeck, K. 1989 Harmonic maps into Lie groups (Classical solutions of the chiral model). 
J. Diff. Geom. 30, 1-50. 

Ward, R. S. 1988 Soliton solutions in an integrable chiral model in 2+1 dimensions. J. Math. 
Phys. 29, 386-389. 

Ward, R. S. 1995 Nontrivial scattering of localized solitons in a (2+l)-dimensional integrable 
system. Phys. Lett. A 208, 203-208. 

Zakrzewski, W. J. 1989 Low dimensional sigma models. Institute of Physics Publishing. 



9 



Figure 2: The soliton configuration at different times for tfie non-generic example f — 1/z 
and h = {z — l)/z^. 
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